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Abstract 

We obtain the effective action of four dimensional quantum gravity, induced by N 
massless matter fields, by integrating the RG flow of the relative effective average ac- 
tion. By considering the leading approximation in the large N limit, where one neglects 
the gravitational contributions with respect to the matter contributions, we show how 
different aspects of quantum gravity, as asymptotic safety, quantum corrections to the 
Newtonian potential and the conformal anomaly induced effective action, are all repre- 
sented by different terms of the effective action when this is expanded in powers of the 
curvature. 



1 Introduction 



Even if we are still lacking a quantum theory of gravitational phenomena we are starting 
to accumulate interesting partial results which will probably be important bits of the final 
theory. These are, among others, the conformal anomaly induced effective action, which was 
first written down in [lj; the low energy corrections to the Newtonian gravitational potential 
analyzed in [2J; the possible ultraviolet (UV) finite completion of the theory described in the 
asymptotic safety scenario [3]. In this paper we want to show how these different aspects can 
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all be seen as arising from different terms in the gravitational effective action when this is 
expanded in powers of the curvature. In particular, the UV properties of the theory are related 
to the zero and first order terms, i.e. to the renormalization of Newton's and cosmological 
constants, and to the renormalization of the coupling constants of all higher order local 
invariants; quantum gravitational corrections to the Newtonian potential are encoded in the 
finite part of the curvature square terms, while the conformal anomaly induced effective 
action is just one of the possible four curvature structures. 

We will use the effective average action formalism to obtain the gravitational effective 
action induced by Nq massless scalar fields, Ni massless Dirac spinors and Ni abelian gauge 
fields, interacting solely with the background geometry, as the result of the integration of the 
RG flow. In the large N expansion!*] one assumes that the number of matter fields N grows 
large. In the leading approximation one simply neglects the gravitational contributions with 
respect to the matter contributions: this eliminates computational and conceptual issues 
related to the treatment of gravitational fluctuations. Within the effective average action 
formalism this point of view has been analyzed in [6J, where the focus was on the local 
part of the effective average action and on the related renormalizability issues. Here we will 
extend the large N analysis to a non-local truncation of effective average action where this 
is expanded in powers of the curvature. 



2 Effective average action for matter fields on curved 
space 

Following [6j we consider massless matter fields on a curved four dimensional manifold 
equipped with a metric g^. The bare action we consider describes iVo massless scalar fields, 
Ni massless Dirac spinors and N\ abelian gauge fields interacting only with the background 
geometry: 



S[(f), ip,An,c, c; g] 



No rl 



Ni 



i=l 



i=l 



Ni M 1 



(1) 



1 N is any of Nq,Ni,Ni 
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Here x is a parameter, only when x = 1 the scalar action is conformally invariant if the 
scalar has conformal weight one. The Dirac operator is defined using the covariant Dirac 
matrices 7^ = e^ a . Note also that on an arbitrary curved manifold the abelian ghosts do 
not decouple and cannot be discarded; we will choose the gauge a = 1 from now on. 

The effective average action Tk[4>, xp, A^, c, c; g] is a scale depended generalization of the 
standard effective action (depending on the infrared (IR) cutoff scale k) that interpolates 
smoothly between the bare action for k — > 00 and the full quantum effective action for k — > 
pjj. It satisfies an exact RG flow equation describing its flow as the IR scale k is shifted. The 
flow equation relevant to the field content we are considering reads as follows: 

ar^rf-itt *^ (2) 

where r£ 2 '°' ) [<£>, g] is the Hessian of the effective average action taken with respect to the 
field multiplet ip = (</>, A^, c, c) and the trace is a "super-trace" on this space. Equation 
(T5]) is exact and is both UV and IR finite. The flow equation ([2]) can be seen as a RG 
improvement of the one-loop effective action which is derived from the modified bare action 
S[(p, g] — >• S[(p, g] + ASk[(p, g], where ASk[<p, g\ is a cutoff action quadratic in the fields. This 
cutoff action is constructed in such a way to suppress the propagation of all field modes smaller 
than the RG scale k. The effective average action and the exact RG flow it satisfies (j2J) offer a 
different approach to quantization: in theory space, the space of "all" action functionals, the 
bare action represents the UV starting point of a RG trajectory which reaches the quantum 
effective action in the IR. The integration of successive modes is done step by step as the 
cutoff scale k is lowered. More details on the construction of the effective average action on 
curved backgrounds or in presence of quantized gravity can be found in [5]. 

As explained in |8j, since we are considering matter fields interacting only with the back- 
ground geometry, we can replace the the Hessian of the effective average action in the rhs 
of (p) with the Hessian of the bare action (Jl]). After performing the field multiplet trace we 
finci: 

dtTkl^^^^^qg] = — Tr h k (A ) - -y-Tri h k (Ai) 



^Tri h k (Ax) - Tr h k (A gh ) 



(3) 



where we defined the function h k {z) = f^ fc ^ and Tr s indicates a trace over spin-s fields. 



2 We are using a Type II cutoff operator in the nomenclature of [5]. 



3 



The differential operators introduced in ([3]) are the following: 

A = A + |i? Ai=A + ^R {A x y v = Ag» v + R» v A gh = A , (4) 

where A = — g^V^Vu is the covariant Laplacian. In the next section we will determine the 
RG flow of the effective average action by calculating the functional traces on the rhs of fl3]). 



3 Flow equations and beta functions 



The traces on the rhs of the flow equation ([3]) can be expanded in powers of the curvature by 
employing the the non-local heat kernel expansion [7J as was done in [HI [9]. To second order 
we find the following result^: 



iV - 4iVi + 2iVi r 
{Alt) 2 d t T k [g\ = ^ Q 2 [h k ] J d 4 xyfg- 

(l-x)N + 2Ni-AN 1 
+- ttt^ Qi[h h ] / d^Xy/gR 
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+ / d 4 x^gR 



ds h k (s) s F R (sA) 



R 



+ / tfxy/gC, 



lMua/3 



ds h k (s) s Fc(sA) 



o 



+ 



N + HJVi + 62iVi 

Qo[h k ] I d d x^E 

N + Ni- 3N! 



60 



Q [h k ] / d d XyfgAR 



(5) 



where 1Z stands for any curvature. The non-local heat kernel structure functions Fq{x) and 
Fr{x) in (jSJ) are linear combinations of the basic non-local heat kernel structure function 
f(x) = J? d^ e~ x ^ x ~t) anc l read as follows: 



5N - 20iVi + ION, f f x ) _ i (3 - 2 X )iV - 2JVi - 2^ ff x ) 
Fr{x) = ^ iA±. + 



48 



48 



x 



;43 - 42 X )iV + 8iVi - 22N, i ( 3 _ 2x ) 2 N - 2N X 



288 



x 



576 



/(*) 



^Hcre we define T k [g] = T k [0,0,0,0, 0; g] 
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N - 4iVi + 2Ni fjx) - 1 _ Ni - 2Ni fjx) 
°^ x ' ~ 4 x 2 Ax 

N + 2N L - lOJVi l jVi 

+ i* ; + f'<*>- < 6 > 

In (jSJ) /ifc(s) is the (inverse) Laplace transform of the function hk(x) and the Q-functionals 
are defined as Q n [f] = f^°dzz n ~ 1 f(z) when n > and as Q n [f] = {— l)/^(0) when 
n < 0. 

Since the RG flow has generated all possible terms compatible with diffeomorphism in- 
variance on the rhs of (jSJ), we need now to consider the most general truncation for the 
gravitational effective average action T k [g] to insert in the lhs of (j3J). As proposed in [HI [9], 
we consider an ansatz where the effective average action is expanded in powers of the cur- 
vature and where the scale dependence is carried both by running couplings and by running 
(possibly non-local) structure functions. To second order in the curvatures the expansion of 
the gravitational effective average reads as follow^: 

1 



F k [g) = j d A x^fg 



(2A k -R) + Rf R , k (A)R+ 



16ttG*. 

' ' o(n 3 ). (7) 



+C, va pfcA&)C^ + -E + -AR 

Pk T k 



In (JTj) the couplings A k and G k are the running cosmological and Newton's constants, fc,k 
and fn jk are the two independent curvature square running structure functions (which at this 
point are arbitrary functions of the covariant Laplacian A) , while p k and r k are the running 
couplings related to the Euler and total derivative invariants. The other two curvature square 
running couplings \ k and £ k are related to the running structure functions by the following 
relations: 

/c,*(0) = r^- /ft*(0) = i (8) 

l; k 

We can extract the beta functions for the couplings A k , G k , p k , r k and the flow equations 
for the running structure functions fc )k , fn, k by comparing (j3J) with ([?]). In particular, by 
matching the coefficients of the operators J yfg and J ^JgR we find the following relations: 

/ I \ (1 - x)N + 2Ni - 4JVi 



4 This implyes also that we are adding to the action (TTJ) the non-dynamical term Fa[<7]. 
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while by matching the non-local curvature square terms we find the following flow equation 
for both running structure functions: 



d t fi,k(x) = / ds h k (s) s Fi(sx) , (10) 

( 47r ) Jo 

where i = C,R and x stands for A. If we insert the explicit form of the functions Fi(x) from 
(JH]) we can rewrite equation (1101) explicitly in terms of Q-functionals: 



{^) 2 d t f Rtk {x) 



(3 - 2 X ) 2 iV - 2N, 
576 

(3 - 2 X )N - 2Ni - 2iVi 

H 

48 

(13 - 12x)iVo + 8iVi - 22^ 



288 

5N - 20Ni + l(Wi f r 1 ) 
f 48^ \J d ^Q^ h ^ z + x ^ 1 -0))-Q2[hk]\ (11) 



FdtQo [h k (z+x£(l -£))] + 
Jo 

d^Qi [h k (z + x^l-m 
Qi[hh] 



(47T) 2 d t f C ,k{x) 



[ d^Qx [h (z+ x e(i-o)] + 

Jo 



4;r 

iV + 2iVi - lOiVi 
4- Q:W 4- 

iV - 4iVi + 2iVi f /-i 



iv - 4/Vi + 2iVi r 1 1 
+ |y ^Q 2 [^(^ + ^(l-0)]-Q2N| • (12) 

Finally, matching the remaining local curvature square terms gives: 

„ /i\ JV„ + lliVi + 62JVi 

<*<>*(*) = *5 

(4vr) 2 9 t (-j = QoN- (13) 

The flow equations ([2]), (ITT]) . (j!2p and (11 3 p completely describe the RG flow of the effective 
average action (J7|). 

The Q-functionals in the flow equations just derived can be evaluated once a particular 
cutoff shape function has been chosen; it is possible to evaluate them analytically if we 
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employ the so called "optimized" cutoff shape function Rk{z) = (k 2 — z)6(k 2 — z). We 
easily find the cutoff shape independent value Qq [hk] = 2 and the results Q n [hk] = ~jjk 2n , 
while the evaluation of the parametric integrals of Q-functional present in (TTT]) and (I12p is 
more involved and we refer to the Appendix of [9J for further details. We can now write 
down the beta functions for the dimensionless cosmological constant, = k 2 Ak, and for the 
dimensionless Newton's constant, Gk = k 2 ~~ d Gu, which follow from ([9]): 



d t A k 
d t G k 



N - ANi + 2N± _ (1 - x)N + 2Ni - AN X _ _ 
-2Afc H — 2 Gk H AfeGfe 



47T 

(1-^0 + 2^-4^ 
2G k + G k 



Gtt 



(14) 



The first terms on the rhs of (fl4"|) represent the canonical scaling of the cosmological and 
Newton's constants, while the other terms represent the running induced by the interaction 
with matter fields. Note also that conformal scalars, i.e. x = 1> do not contribute to the 
running of Newton's constant. Equations (TTTj) and (|T2|) can be compactly rewritten as: 



d t fi,k(x) 



(4tt) 



x 

;<)> \ ]7i 



(15) 



where we defined the following cutoff shape dependent functions: 



9i{u) 



ds hkis) s Fj(sk 2 u) . 



(16) 



In (TLB"]) u = x/k 2 stands for the covariant Laplacian in units of k; note that the UV regime 
corresponds to small values of u while the IR regime corresponds to large values of u. The 
integrals in ( TL6l) can be evaluated analytically |9j, we find the following forms: 



9c{u) 



9r(u) 



N + 6A^i + 12iVi 

2 



120 



A^o + 6Ni + 12iVi 8Aq + 8Ni - 64 A^ 

2 2 



120 



120m 



16Aq - 64A^i + 32Ni 



120m 2 



'l - -d(u-A) 



u 



[1-X) 2 N 



[l- x yN (l-x)N + Ni-2N 1 



72 I 72 

N - ANi + 2N 

H 2 - 

18u 2 



18m 



'l - -6(u-i) 
u 



(17) 



Note that the form of ( TT71) implies that conformally invariant matter, contributes to gii(u) 
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and thus to the flow of fn,k{x), showing that for k 7^ the flow generates non-conformally 
invariant interactions. We will see later that the conformal invariance of the curvature square 
terms will be partially restored in the IR at k = 0. The constant terms in (TTT1) . when matched 
with (151) , give the beta functions for the couplings and These, together with the explicit 
forms of the beta functions f lT5|) . are: 

1 N + 6N1 + 12iVi 
dAk ~ "(4^ 60 Afe 

1 (1- X ) 2 iy 2 

1 jV + HiVi + 62JVi 2 
^ = (4^ 180 Pfc 

1 iVo + Ni - 3A^ 

S ' T ' = "(S? i * < 18 > 

The numerical coefficients in (Tl8|) are scheme independent, i.e. they don't depend on Rk{z) 
and on the other details of the cutoff choice. Note that the beta function of the coupling 
vanishes when one is considering conformally invariant matter. The beta functions ( [141) and 
(TL8]) agree with those obtained in j6]. 



4 Asymptotic safety 

We proceed now to integrate the RG flow of the couplings of the local curvature invariants 
present in the truncation we are considering, i.e. the cosmological, Newton's, and the cou- 
pling constants of the curvature square terms. We will see that precisely the flow of these 
couplings is related to the UV renormalization of the theory, that in the effective average 
action formalism is encoded into the initial conditions of RG flow. These relations will tell 
us how to choose the bare couplings in order to achieve a finite continuum limit when the 
UV regularization is removed. 

Since (Tl4|) is a closed system for and Gk, we can start to look for non-Gaussian fixed 
points of the RG flow, that we need to find in order to construct a continuum limit, by first 
solving: 

d t K k = d t G k = 0. (19) 

It's easy to see that the system f[T9"|) admits both a Gaussian fixed point = Gk = and a 
non- Gaussian one [B]: 



S 



Figure 1: Renormalization group flow in the (A^G^) plane form equations (124")) and 02 8p . 



A = 3 N -Wt+2 Nl = Un 

4:(l-x)N + 2Ni-4:N 1 * (l-x)N + 2Ni-AN l ' 1 ; 

Note that to have a fixed point with positive Newton's constant we need to satisfy the 
inequality (1 — x)Nq + 2Ni — 4Ni > 0. The stability matrix around the non-Gaussian fixed 
point (1201) is readily calculated and shows that it is attractive in the UV 0. This shows that 
both the cosmological and Newton's constants are asymptotically safe couplings. This fact 
shows that, within the the truncation we are considering and in the large N limit, four 
dimensional quantum gravity is asymptotically safe§. The beta function system (Il4p can be 
easily integrated analytically. Integrating the second equation in (I14p from the UV scale A 
to the IR scale k gives the following relation connecting the fc-dependent, fixed point and 
bare Newton's constant^: 

1 1 / 1 l\/A\2/7„\ 2 

(21) 



1 


1 






®" 


fk 


Gk 


G* 


{g a 


GJ 







In (|21|) we introduced the reference scale ko, which will play a role similar to the renormal- 
ization scale of standard perturbation theory. To make (l2"Tj) finite in the continuum limit, 



5 It has eigenvalues —4,-2. 

6 See [5] for a discussion of asymptotic safety when also gravitational fluctuations are considered. 

7 If we define the couplings = and B^ = 16 ^ Gk the system (fT4")) becomes d t A k = 4(A* — Ak) 

and dtBk = 2(B* — Bk) where Ak — Akk 4 and Bk — Bkk 2 . These relations immediately give Au = 
A* + (Aa - A*) (x) 4 and B k = B* + (B A - B*) (f ) 2 from which we read off dHJ) and dHJ). 
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\G A 


GJ 





i.e. in the limit A — > oo, we need to renormalize Newton's constant; this can be done by 
imposing the following condition: 

' 1 1 \ / A x 2 

C G , (22) 

where Cq is a constant that we will determine later. Solving (122]) with respect to bare 
Newton's constant G A gives: 

G 

G A = -* T. X 2 > (23) 

1 + CgG* (^) 

which shows how we need to choose G\ to approach G*, as A — > oo, in order that the lhs of 
(122]) remains constant. Inserting now (122]) in (12T]) gives the functional form of the ^-dependent 
Newton's constant: 

G* 

1 + CgG* (%) 

Equations (123]) and (liH]) are formally the same at the level of the approximation we are 
considering but have different physical meanings: the first equation tells us how the bare 
coupling changes as we vary the UV scale A, while the second equation shows how the k- 
dependent coupling flows as we integrate more and more degrees of freedom by lowering the 
RG scale k. From (124j) we see that the dimensionful Newton's constant, Gk = k~ 2 Gk, reaches 
the following renormalized value at k = 0: 



Gk = ^T7T2 ■ (24) 

') 



G ° ~ 7^T2 =^ °G~ 77T2 > (25) 



which fixes the value of the constant introduced (|22|) . Similarly, we obtain the following 
relation between the /c-dependent, fixed point and bare cosmological constants: 

After we impose the UV renormalization condition: 



' Aa _ A,\ ( A- 
G A GJ Uo 



i 



C A 
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which can be solved for the bare cosmological constant: 



Aa 



A, + CaG* (g) 



(27) 



we obtain the ^-dependent cosmological constant: 



A fc 



k+cag, (|y 

1 + 



(28) 



As before, f l27]) and f )28|) are formally equivalent. The renormalized value of the dimensionful 
cosmological constant = fc 2 A& is obtained from (J2SJ) and shows that Ca = ^rp- Since 
both Gfc and are relevant parameters in the UV (they are both attracted by the non- 
Gaussian fixed point) we need to fix their value from experiments, i.e. we need to measure 
Ao/fcg an d Go&o- The phase diagram obtained by plotting equations ( 124)) and ([28]) is shown 
in Fig. 1; we can clearly see that all flow trajectories reach the IR without any obstruction, 
a feature that is not yet present in the analogous RG flows that include also gravitational 
fluctuations [5]. Note that the dimensionless product A Go is independent of the arbitrary 
scale fco, thus AoGo is a true IR numerical prediction associated to a given RG trajectory. 
The Einstein-Hilbert terms in our ansatz (J7]), in the k — > limit becomes: 



Obviously this part of the gravitational effective action is not conformally invariant due to 
the presence of the scales A and Gq. We see from ff2T]) and (126]) that the only way to obtain 
conformal invariance in the IR is to set the UV couplings to their fixed point values: A* = Aa 
and Ga = G* so that Ao = and 1/Gq = 0, obviously the only way to restore conformal 
invariance at the quantum level. 

Next we can integrate the beta functions (fl8j) for the couplings \k,£,k, Pk, T k- These cou- 
plings have only the Gaussian fixed point, at least within the approximation we are consid- 
ering [10] . This is easily done: 



where gk stands for any of the couplings \k,£,k, Pk, Tk and where the relative constants C can 




(29) 




9k 9a k 



— + C log A - C log ^~ 



9k k k 
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be read-off from ( I18jl . We can now choose the following UV boundary conditions: 



1 



(30) 



9 a = 




which imply the following ^-dependence of the coupling constants: 



1 



(31) 



9k = 




As we mentioned before, (13"U|) and (T3"T|) are formally equivalent; they both show that these 
couplings are asymptotically free. Note that when we reinsert the gk from (}3~Tj) in the effective 
average action ([7]) and we take the limit k — > the local curvature square contributions 
diverge. It seems that the effective action has an IR divergence. We will see in the next 
section that this is not the case: when we include the finite curvature square contributions 
in ([7]), the logk in (}3~Tj) combines with another such term, coming from the running structure 
functions, so that their sum has a finite k — > limit. 

5 Curvature square terms and corrections to the Newto- 
nian potential 

In this section we start to look at the finite non-local part of the effective action and we show 
how this part of the action can be obtained as the result of the integration of the RG flow. 
We focus on the running structure functions in the curvature square terms, where we can 
perform all steps analytically. 

By integrating the flow equation ([15]) for the running structure functions from the UV 
scale A to the IR scale k we find the relation: 



where fi.h_{x) are the bare structure functions, while fi t k(x) are the /c-dependent ones. Chang- 
ing variables to u = x/k 2 , with dk/k = —du/2u, gives: 



It should be noted that when we insert the explicit forms of the functions gi(u) from ffTTl) . 





(32) 
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the constant terms will make the integrals in (152]) logarithmically divergent at the lower limit 
when A — > oo. We can isolate these divergencies by subtracting these constants from the 
integrals in (152]) in the following way: 



(4vr 



fc,k{ x ) = fcA x ) 



1 N + 6iVi + 12iVi k 
— log 



1 iV + 6N1 + 12JVi A 

log — 

(4tt) 2 120 S k 



1 



x/k2 du 



N + 6N L + 12Nt 



9c (u) 



and 



120 '"° k (4tt) 2 J x/A 2 2u 

fR,k{x) = f R ,A{x) ~ j—y 2 



120 



(33) 



log 



+ 



1 (1-X) 2 N k_ J_ 
(4vr) 2 72 ° g k (4vr) 2 



x/k2 du 



x/A 



2 2u 



9R W 



A 

(1 - x) 2 N 

72 



(34) 



The log y~ terms in the first lines of (1331) and (|34p are the UV divergences that are also 



encountered in standard perturbation theory [TT]. We can renormalize the running structure 
functions by imposing the following UV boundary conditions: 



1 A^o + 6Ni_ + 12JVi A 

fcAx) = lio hg Y Q +Cc 

1 (1 - X ) 2 N , A 
fRA{x) = (47r)2 72 hg k- +CR > 



(35) 



that make the first lines in (l33l and (1341) vanish. Here cc and Cr are possible finite renormal- 
izations constants. In this way the second lines of (133]) and f }34"|) are finite in the UV A — > oo 
and in the IR k — > limits. Note also that the first terms of the second lines of ( 133]) and ( |34j) , 
when combined with (jSJ) correctly reproduce (13T]) . Performing the integrals in (133]) and ( 1541) 
now gives the following finite contributions to the curvature square part of the renormalized 
effective action: 



To foil 



1 



7v 2 



2(4tt) 



/ 



d 4 Xy/g 



N + 6Ni + 12JVi 



120 



-C. 



log ^C' i ^ /3 

K 



+ (1 -5/ JV ° j i iog^ 

/Cn 



72 



(36) 
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where we fixed the renormalization constants to the values: 



23 Ao 
1800 



3JVi 



i AN, 



1 + 10(1 - X ) + 30(1 -x) 
2160 



2 A ( 



o 



Ni 



cc = 



cr = 



+ 



2 



50 



75 



360 



120 



The effective action (136]) is the same as the one obtained using standard perturbation theory 
if we equate the arbitrary scale k$ with the renormalization scale /x [H]. In ( 136]) we dropped 
the Euler and the total derivative terms, as we mentioned in the previous section, these terms 
are not finite in the IR limit k — > since their running is of the form fl3T]) . but this does not 
cause any problem since these terms are total derivatives and thus vanish in genera ^. Note 
also that if we were employing a different cutoff shape function Rk{z) than the one employed 
here, we would had found a different form for the flow for k ^ but, as was shown in the 
two dimensional case [8], the k — > limit would had always been equal to (136]) . 

In ( 136]) the only contribution to the renormalized Ricci structure function fRfi(x) comes 
from the scalar fields and this contributions vanishes if we consider conformally invariant 
scalars by setting 1. Even if x = 1 the functional (136]) is not conformally invariant. This 
happens because we are considering an ansatz of the effective average action structured as in 
([7]), which by construction cannot accommodate a conformally invariant limit for k = since 
the expansion we are performing is not so. As we will explain in the following section, we 
expect that when % = 1 the only non-conformally invariant part of the effective action to be 
the conformal anomaly induced effective action |14j . The form found in equation ( 136]) must 
thus be the result of the expansion of a more general conformally invariant term, that to 
order curvature square reduces to (136]) . We can speculate such a term to be of the following 
form: 



where the operator Oc is a fourth order conformally invariant differential operator that acts 
on four indices tensors with the symmetries of the Weyl tensor. The explicit form of this 
differential operator is unknown [12] but it must have a form such that (1371) has a smooth 
flat space limit. At this point it is important to note that a conformal variation of ( 136]) . when 
X = 1, will generate the anomaly term proportional to C 2 with the correct coefficient (see 
equation ( l4~4"]) ). This fact has already been noticed by several authors [TT1 [13], in particular 
the second authors consider this a genuine anomaly term in contrast to the expectation that 
effective action proposed in [lj is the only one that produces the conformal anomaly upon 

8 The integral f y/gE = 32ir 2 x is the Euler characteristic of the four dimensional manifold and is thus 
non-zero in general. The renormalization of this term is probably related to topological fluctuations. 




(37) 
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variation [M]. Since these first authors obtained the form (1361) by a flat space calculation, we 
can think of them to be in the same position as we are in our calculation: the non-conformally 
invariant form, and thus anomalous term, is what remains of a truly conformally invariant 
term of the form (l3~Tj) that reduces to a non-invariant term when, to be able to perform 
calculations, an expansion that brakes conformal symmetry is performed. This point of view 
can compose the different expectation expressed on one side by [13] and on the other side by 
1 14] and thus deserves further study. 

A physical reason to expect a term like (1361) in the effective action, at second order 
in the curvatures, is that it generates the quantum gravitational corrections to Newton's 
gravitational potentials These corrections can be obtained along the lines of [9j and are the 
following: 



The first part of the quantum gravitational correction in ([3*5]) is the purely gravitational 
one found in [9]. As explained in this last reference, the correction (I38p is equivalent to 
the one obtained using effective field theory arguments. Note that every field species (even 
conformally invariant matter) gives a positive contribution, of the same form of the purely 
gravitational one. From ( )38|) we can see that when the number of massless matter fields N 
grows large, their contributions soon dominate the purely gravitational contribution. 

6 Higher order terms and the conformally induced effec- 
tive action 

In the last two sections we explicitly derived and integrated the RG flow of the gravitational 
effective average action to second order in the curvature expansion. In this section we now 
look at the third and fourth order terms. 

One can consider the curvature cube terms in the effective average action ansatz (J7J) and 
can project their flow by employing in (~~~) the non-local heat kernel expansion to third order 
[15J. There are now ten independent running structure functions and the curvature expansion 



3 This is only the "graviton polarization" part of the correction. 
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of the effective average action has the following form: 



i iU r 

rk[g}\w = (^pE J d d x^/gf i!k (A h A 2 ,A 3 )KiK 2 n 3 (i) , (39) 

where the curvature structures lZxlZ 2 lZ 3 (i) are given in [15]. For example 1Zi1Z 2 1Z 3 (1) = 
R1R2R3, Hi 7^2^3(2) = Rf^R 2u R 3a and so o i0; following the conventions of [15] the Lapla- 
cians Aj act only on the curvatures TZi for i = 1,2,3. The third order terms on the rhs of 
the flow equation (j3J) can now be written as: 

j 10 poo _ r 

d t V M\n* = J^ylYlJ Q dsh k (s)s^J d d x^gF i (sA 1 ,sA 2 ,sA 3 )Tl 1 Tl 2 Tl 3 ({), (40) 

where the explicit form of the non-local heat kernel structure functions Fi{x\, x 2 , x 3 ) can be 
obtained from those given in [15] . Since the curvature structures lZilZ 2 lZ 3 (i) may contain 
covariant derivatives, third order terms are proportional to different powers of s. In fact we 
have a\ = a 2 = a 3 = 3, a 4 = a 5 = a e = aj = 4, a% = a 9 = 5 and aio = 6. By comparing 
fl39|) to (I40p we obtain the flow equations for the third order running structure functions: 

d t fi, k {x u x 2 ,x 3 ) = j^ gi {—,-,^ , (41) 

where we introduced the scheme dependent functions: 

1 f°° 

g i (u 1 ,u 2 ,u 3 ) = dsh k (s)s a '~ 1 F i (sk 2 u 1 ,sk 2 u 2 ,sk 2 u 3 ) , (42) 



(An) 



where Ui = Xi/k 2 for i = 1,2,3. At this order, computations become increasingly hard and 
it seems rather difficult to straightforwardly evaluate the s integrals on the rhs of ( l42l) . even 
is we employ the optimized cutoff shape function R k (z) = [k 2 — z)9(k 2 — z). It was shown 
in [6] that if one employes this cutoff shape function, the couplings of the cubic and higher 
invariants, i.e. gi tk = ^^(0,0,0), have zero beta functions. With this cutoff choice, we need 
to UV renormalize only the couplings we already considered at lower order of the curvature 
expansion. In spite of this, the structure functions gi(xi, x 2 , x 3 ) are non-zero even if they have 
a vanishing Taylor expansion around the origin (since Q- n [h k ] = for all n >1), probably 
they are proportional to theta functions as are the second order functions gi(u) (|T7|) . For 
this reason no UV divergencies arise when we integrate the flow equation (T42|) and we find a 



in 



Note that we are considering only the purely gravitational curvature structures of |15| . 
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finite result for the renormalized third order structure functions: 

1 f°° dk 



fifi(xi,X2,X a ) 



(4tt) 



k 



■9i 



fc 2 ' /c 2 ' fc 2 



(43) 



In (l4*3|) we imposed the UV boundary conditions fi t \(x\, X2, £3) = 0. The absence of UV 
divergencies can also be understood by looking at the local expansion of the gravitational 
effective average action where, as we just noticed, the dimensionful beta functions turn out 
to be zero d t gi^ = 0, thus implying = (^j 1 g^ where g^ = g^ are the dimensionless 
coupling with dj < 0. Thus in the limit A — > oo is finite and we have g^ — > 0. We finally 
obtain the cubic curvature part of the gravitational effective action r [g]|.£3 when we re-insert 
flU into PJ). 

It will be interesting to extend the curvature expansion to the next (quartic) order in the 
curvatures, since the conformal anomaly induced effective action is of this order pQ. This 
effective action is usually obtained by integrating the conformal anomaly; for the field content 
we are considering and for x — 1 this is the following: 



(47T) 1 



iVo + 6iVi + 12JVi N + lliVi + 62 N x 

1 C 2 rr 1 E 



120 

N + Ni 



360 



3N X 



+ 



60 



-AR 



(44) 



To derive 



one uses the relation (Tff 



(47r)" 2 / v /^tr6 4 (S ,(2) ), where S i2) is the Hessian 
of (00) and the 64 are the relative local heat kernel coefficients. It was shown in [lj that the 
following effective action has an energy-momentum tensor whose trace reproduces ( 1441) : 



^CA[g] 



8(4tt) 



d A x^/g 



N + 6N1 + 12JVi 



N, 



t- HJVi + 62A^ 

2 



360 



60 



E + -AR) 



-C 2 



A 



4,0 



E 



-AR) 



(45) 



where A 4 = A 2 + 2i? / " / V M V 1/ + f-RA + iV^-RV^ is the conformal covariant fourth-order 
operator acting on conformal invariant scalars. The conformal anomaly induced effective 
action ( |44|) will constitute, together with other conformally invariant terms, the four cur- 
vature terms in the curvature expansion of the effective action. But since to quartic order 
in the curvatures we don't even know the form of the non-local heat kernel expansion, any 
calculation along the lines of the previous sections is, for the moment, a very hard task to 
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perform. In spite of this we believe that the integration of the flow of the effective average 
action will generate the non-local effective action ( )44l) as one of the many fourth order terms. 
In [8] it was shown that this is indeed what happens in the analogous two dimensional case. 

7 Conclusions 

As a first step toward understanding the gravitational effective action we studied the effective 
action induced by massless matter fields interacting solely with the background geometry. We 
proposed an expansion of the effective action in powers of the curvatures and we showed how 
various known results about quantum gravity could be seen as arising from different terms 
of it. We showed that asymptotic safety and UV renormalization are related to the running 
of the couplings of local curvature invariants, the first of these being the cosmological and 
Newton's constants. Then we derived the quantum gravitational corrections to the Newto- 
nian potential from the quadratic curvature terms and we explained that all other quantum 
gravitational effects, including those induced by the conformal anomaly are contained in the 
higher order contributions. 

Once matter contributions are well understood, one can infer something about the grav- 
itational effective action when matter contributions become dominant in the large N limit. 
When instead gravitational fluctuations become large, we must consider the fact that gravity 
may actively reacts to matter inducing corrections to the effective action that may be of 
the same order as those induced by matter. The exact solution of two dimensional quan- 
tum gravity shows that this is indeed what happens in the two dimensional case [16]. It is 
thus important to extend this work and the program started in |9j to the treatment of full 
gravitational dynamics. 
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